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ABSTRACT 

In this proof-of-concept study we demonstrate that in a binary system mass can be transferred 
toward an accreting compact object at extremely high rate. If the transferred mass is efficiently 
converted to X-ray luminosity (with disregard of the classical Eddington limit) or if the X-rays are 
focused into a narrow beam then binaries can form extreme ULX sources with the X-ray luminosity 
of Ax ^ lO^^ergs”^. For example, Lasota & King argued that the brightest known ULX (HLX-1) 
is a regular binary system with a rather low-mass compact object (a stellar-origin black hole or a 
neutron star). The predicted formation efficiencies and lifetimes of binaries with the very high mass 
transfer rates are large enough to explain all observed systems with extreme X-ray luminosities. These 
systems are not only limited to binaries with stellar-origin black hole accretors. Noteworthy, we have 
also identified such objects with neutron stars. Typically, a 10 Mq black hole is fed by a massive 
10 M©) Hertzsprung gap donor with Roche lobe overflow rate of ~ 10“^ Mq yr“^ (« 2600MEdd)- 
For neutron star systems the typical donors are evolved low-mass 2 M 0 ) helium stars with Roche 
lobe overflow rate of ~ 10“^ Mq yr“^. Our study does not prove that any particular extreme ULX is 
a regular binary system, but it demonstrates that any ULX, including the most luminous ones, may 
potentially be a short-lived phase in the life of a binary star. 

Subject headings: stars: black holes, neutron stars, X-ray binaries 


1. INTRODUCTION 

Our Universe is populated with black holes (BH) and 
neutron stars (NS) in various binary configurations. In 
our Galaxy, many such binaries show significant X-ray 
activity suggestive of a mass transfer and accretion onto 
the compact star. X-ray luminosities of only two Galac¬ 
tic X-ray binaries (XRBs) exceed ~10^®ergs“^ (approx¬ 
imately th e Eddington limit for a lOMfn BH) in GRS 
I915-I-105 (|Fender & Belloni||2004|) and possibly i n SS 
433 if the system were observed along the jet axis (Fab- 
rika et aL]|2006 1. 


However, a large population of extra-galactic point-like 
X-ray sources with lumin osities in excess of erg s~^ 

has been identified (e.g., Fabbiano et al. 1989 Liu||2011 


Walton et al. 

2011) 

ray sources ( 

JLA.j 


These so-called U 


1989 _ 

tralummous X- 


tion onto a supermassive BH (M > 10® Mq) can be ex¬ 
cluded as the source of their luminosity. Instead, the 
two most popular scenarios to explain their nature in¬ 
clude binary systems hosting (i) a, stellar mass BH ac¬ 
creting at a super-Eddington rate, or (ii) an inter mediate 
mass BH (IMBH) a nd sub-Eddington accretion (Golbert 
& Mushotzky|1999 ). In the latter case, formation ot HHs 
heavier than ~ 100 Mq presents a problem for current 
models of stellar evolution. Although, it has been sug¬ 
gested th at such IMBHs may be fo rmed in dense globular 
clusters (Miller & Hamilton 2002) or even as a result of 
stellar evolution ot very massive stars ( ^ 200 — 300 Mq; 
Growther et al.j 2010 Yusof et al.||2013). 

'i'he super-Eadington BH accretion invoked in the stel¬ 
lar origin scenario remains still a relatively poorly under¬ 
stood regime. However, several theoretical mechanisms 


able to breach the Eddington limit have been proposed, 
e.g., beaming and/or hyper-accretion allowed by non- 


bubbles”) (King et al. 

2001 

Begelman et al. 12006 

Pouta- 

nen et al.||2007 

King 

2008 

, as well as a contribution of 


ULX population ( Medvedev fc Poutanen|[20l3 ). 


Robust observational constrains on the mass of the 


accretor in some ULXs (e.g., Motch et al. 2014) indi¬ 
cate that indeed super-Eddington accretion ont o stellar- 
origin BHs i s realised in nature. Moreover, [Bachetti] 


et al. (2014) investigating the M82 X-2 source have 
demonstrated that the super-Eddington accretion is also 
possible in XRBs hosting a NS. Additionally, some 
sources transit relatively fast between the super- and sub- 
Eddington regi mes on timescales as short as a few day s 
to a week (e.g., Walton et al.|2013 Bachetti et al.|20f3 ). 
Such shor t timescales are in contradiction with IMBH 
accretors ( LasotapOlS ). 

Nevertheless, it has been speculated that the bright¬ 
est ULXs, with lumin osities > 10"^® ergs ~^ may be 
candidate IMBHs (e.g., |Sutton et al.||20i2|). Recently, 
compelling evidence, based on quasi-periodic oscillations 
(QPO), was presented in sup port of a ~ 400 Mq BH in 
M82 X-1 (Pasham et al.|2014). However, it was also sug¬ 
gested that these QPUs may be harm onics of puls ar rota¬ 


tion period (Kluzniak & Lasota|2015 ). Moreover, Sutton 
et al. (2015) demonstrated thaf the IMBH candidate m 
IG 4320 is actually a backgro und AGN. To date HLX-1 
with Lx = 1.1 X 10^^ ergs“^ ( Farrell et al.||20fT ) is the 


brightest known ULX (for a d iscussion ot the brightest 
ULXs see Servillat et al. 2011). 


We approach the ULX issue from the standpoint of one 
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particular evolutionary model for binary evolution. We 
consider only the far end of the ULX luminosity space, 
Lx ^ lO^^ergs”^, and we refer to the sources poten¬ 
tially able to reach these luminosities as extreme ULXs 
(EULXs). We explore the possibility that EULXs are 
binary systems with Roche lobe overflow mass transfer 
rates that highly exceed classical Eddington limit. Eor 
the purpose of this proof-of-concept study, we assume 
that the transferred mass is efficiently accreted onto a 
compact object and converted to X-ray luminosity in the 
full range of possible mass accretion rates. This is in 
contrast with the generally accepted view that the con¬ 
version efficien cy decreases with incr easing mass accre¬ 
tion rate (e.g., Poutanen et al. 2007|. However, if mass 
is lost during accretion process and even if conversion 
into X-ray luminosity is not fully efficient, a geometrical 
beaming can provide large X-ray luminosities for sources 
considered in our study. 

2. MODEL 

We empl oy the binary populatio n synthesis code, 
StarTrack dBelczynski et al~ 2008a), with updates as 
discussed in Uominik et al.| (|2T)M ). \Ve em ployed follow¬ 
ing initial conditions: Kroupa et al. (1993) broken power 
law for initial mass function (primary m ass between 
6 - 150Mf7^ ), flat mass ratio distribution (Kobulnicky 


et al. 2006 secondary mass between 0.08 - 150 M~ 


hat in logarithm distribution of separations (Abt 1983 
fja) 1/a) and thermal distribution of eccentricities 
( Duquennoy fc Mayor||199l| /(e) ~ eV 

We use the BDINC platforiiQ for volunteer 
computing in our program ‘‘UniverseQhome ’ ’ 
(http://universeathome.pl) to obtain a large num¬ 
ber of X-ray binaries (N = 10® of massive binary 
systems were evolved by volunteers). The X-ray binary 
is defined as a system hosting a donor star transferring 
mass via Roche lobe overflow (RLOF) to a compact 
object companion (NS or BH). For any given system 
our evolutionary models provide the donor RLOF mass 
transfer rate, Mrlof- The accretion rate onto the 
compact object, Mace, is estimated in three different 
ways from the mass transfer rate (see below). 

We search for evolutionary channels that allow for the 
formation of binaries with highest possible RLOF mass 
transfer rates onto stellar-origin BHs and NSs. Conser- 
vatively, we model o nly solar metallicity {Zq = 0.02; 


Villante et al. 2014) and we allow initial mass func- 
tion (IMl') to extend only to ISOM©. Note that lower 
metallicity and higher mass stars may for m massive 


(> lOOMrri) black holes in binary systems dBelczynski 
etalllMdt . 


To estimate RLOF mass transfer rate Mrlof we eyalu- 
ate donor properties and yarious angular mo mentum loss 
mechanisms i n a giyen binary as described in |Belczynski| 
et al.| (2008a). It is highly uncertain how RLOF mass 
transier rate is to be conyerted to X-ray luminosity, and 
we discuss this issue below. 

The grayitational energy of the RLOFing material 
becomes converted into radiation in an accretion disk 


forme d a round the com pact object (Shakura & Sunyaev 


1973 see Lasota (2015) for a recent review on accretion 
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Fig. 1. — Most important processes behind conversion of RLOF 
mass transfer rate Mrlof info X-ray luminosity. Our parametri- 
sation of these processes is discussed in detail in ^ and used in 
Eq.[T| 


disk physics). There are a number of effects that have a 
pivotal influence on the Mrlof to Lx conversion process 
in high mass-accretion rate disks (see Fig. [^. 

The first conundrum to consider is the role of winds 
launched from the disk surface. Such winds have been 
ubiquitously detected in high mass accretion rate disk - 
dominated states of XRBs (e.g., Ponti et al.|2012 2015). 


Winds are able to remove a substantial fraction of matter 
that is being transferred through the disk toward the 
compact object. At the same time, winds carry away the 
angular momentum and therefore influence the orbital 
evolution. We introduce the quantity /acc to describe 
the fraction of Mrlof that is not affected by the disk 
winds (Mace)- 

Secondly, not all of the photons produced in the vicin¬ 
ity of the accretor get emitted from the disk. Some of 
them are dragged by the inflowing matter and fall onto 


effect (e.g., Ohsuga|2007a hereafter 007 Abramowicz & 

Straub||20081 

iNarayan et al.| 20121 |Sqdowski et al. 12015 

hereafter 
To indue 
that may 
Finally 

515 

ie H 

low 

we 

), the observed X-ray luminosity is reduced. 

lis and other currently unknown processes 
er Lx we introduced the f~ parameter, 
took into account the processes that may 


increase the observed Lx due to non-isotropic emission. 
We utilise the /^ f actor which in addition to the beaming 


^ http://boinc.berkeley.edu/ 


(King et al. 2001) may also include a contribution due 
to disk geometry, black hole spin, column accretion in 
magnetised NS, etc. The beamed emission will always 
exceed the isotropic one (/“*■ > 1; as an example, if the 
emission goes into a cone of opening angle 1°, 10° or 
100°, it will correspond to f~^ ~ 2.6 x 10^, ~ 260, and 
^ 2.8, respectively). However, a beamed source will be 
visible only from the directions enclosed by the cone, and 
thus its detection probability will be lower than that of 
an isotropic emitter. 

The XRB X-ray luminosity is calculated as 
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Lx = f-f 


_ , eGMacc /acc-^RLOF 


R. 


= V^rlofc^, ( 1 ) 


where the radius of the accretor, i?acc, is 10 km for a 
NS and 3 Schwarzschild radii for a BH, e gives a conver¬ 
sion efficiency of gravitational binding energy to radia¬ 
tion associated with accretion onto a NS (surface accre¬ 
tion e = 1.0) and onto a BH (disk accretion e = 0.5), 
rj = /accbo is the total efficiency of the accretion 

flow, and rjQ is the radiative efficiency of a standard thin 
disk. 0 

In our simulations we define a potential EULX source 
as the one with Lx > 10^^ erg s“^. 

2.1. The reference model 

In our reference model we assume: no winds {face = l)j 
no photon trapping (/“ = 1), no beaming (/■*" = 1), and 
always efficient mass accretion rate-into-luminosity con¬ 
version [rjQ as in the standard thin disk). We adopt this 
condition even for very high RLOF mass transfer rates. 
This is quite arbitrary, however, our goal was to estimate 
the highest luminosities potentially reachable by an X- 
ray binary system without invoking strong beaming. On 
the other hand, r] = rjQ may also correspond to a situa¬ 
tion when a significant part of the Mrlof is lost from 
the system {face < 1) or Lx is lowered due to photon 
trapping (/“ < 1) but the beaming compensates these 
effec ts (/+ = l//acc/~)- Numerical simulations (e.g., 
S15 ) suggest that face is indeed small. However, /+may 
reach 10^-10"^ for the cone opening angles of the order of 
a few degrees. 

Therefore, with our reference model we obtain the 
maximum potential X-ray luminosity of a given binary 
system if the non-standard accretion disk effects are neg¬ 
ligible or compensate each other. 

2.2. The Ohsuga model 

Next, we have considered a BH accretion case with 
Mace < Mrlof (.fa ce < 1), and Mace Constrained fol¬ 
lowing the results of |O07| who performed global, axisym- 
metric simulations of super-critical disks with the aP 
viscosity, including effects of outflowing winds and pho- 
ton trapping. We use the parametrisation of the 
results {Mace as a fu nction of Mrlof) derived in 
czynski et al. (2008b), 


007 


Bel- 


Ohsuga (2007b), the corresponding /+ will be equal 66. 
The accretion onto NSs is limited to th e classical Edding¬ 
ton limit in this model ( Ohsuga|2007b ). Lx is calculated 
as in Eq. We note that we have extrapola ted results of 
the original calculations of Ohsuga (2007a|) to the range 


of mass transfer rates we have in our simulations. 

2.3. The Sqdowski model 

Finally, we have experimented with the cons train ts on 
the Mace and Lx following the recent results of S15 The 
accretion rate at the BH horizon was found to be only a 
fraction of the RLOF rate. 


Mace = Mrlof 


Rin 

RciU 


(3) 


where Rout and Rin = 20Rg (Sqdowski et al. 2014[) 


are the outer and inner radii of the wind emitting re¬ 
gion, respectively. It is hard to estimate the location 
of the outer edge, so we assumed a constant fraction 
Rwind/Rout = 0.01 (the wind probably is not emitted out 
to the edge of the disk), which corresponds to face = 0.01. 
In this model we utilised a different formula for X-ray lu¬ 
minosity which already accounts for the effects of photon 
trapping and beaming, 


Lx = 4 X 10^® X e 


Ma„ 


MEdd M© 


Mbh _i 
ergs , 


(4) 


where 9 is the viewing angle, MEdd = 2.44 x 

10^®;j^gs“^ is the Eddington accretion rate. In our 

simulations we incorporated 9 in the range 0° - 30° (the 
opening angle is 60°). Th is will correspond to f ~*~ « 8 
in our reference model. In Sqdowski et al. (2015) model 
we get different luminosities tor different viewing angles, 
which was included in our simulations to calculate prob¬ 
abilities of obser ving a particular system as the EULX. 
Even though the S15 model was constructed for systems 
with BHs, we assumed that the same prescription is valid 
for the NSs accretors. We note that the results of 10071 
can be put in the framework of this model by taking 
-Rout ~ 50Rg, which is approximately the effective circu¬ 
larization radius (and the outer edge of the disk) of the 
gas injected into their simulation box. We no te that the 
Eq. |4|was obtained by|Sqdowski et al.| (|2015|) for super- 
massive BH, but we extrapolated it to the stellar mass 
ones. 


log 


Ma.. 


Me, 


= 0.934 log 


I Mrlof 




- 0.380. (2) 


crit 


where Merit = 2.6 x 10 ®(Mbh/10 Mq) Mq yr ^ is the 
critical mass accretion rate. To obtain X-ray luminosi¬ 
ties we utilised the Eq. using Macc as provided by 
Eq. 0 and substituting /^/+/accM rlof = Mace- For 
a IOMq BH accretor we obtain 2.8 times lower lumi¬ 
nosity if Mrlof = 10“^M yr“^ and 6.7 times lower if 
Mrlof = O.IMyr“^. For illustration, if the collimation 
angle of the outflow is ~ 20°, as in NS simulations of 


^ In our case r]Q gjj = 1/12 and riO,NS 
|Sunyaev|1973|. 


0.2 (e.g., 


Shakura & 


3. RESULTS 

Under the least restricting assumption (no limit on ac¬ 
cretion rate, ly = rjo) our EULX rate/number estimates 
are to be considered upper limits. We find that 1 per 44 
billion binaries could potentially be an EULX with a BH 
accretor, and 1 per 44 billion binaries could potentially 
be an EULX wit h a NS accretor. This estim ate employs 
a canonical IMF (Kroupa & Weidner||2003 ). 

Figure!^ shows our upper limits tor the number of EU- 
LXs withBH and NS accretors. Note that in our ref¬ 
erence model we obtain virtually as many potential EU- 
LXs with a NS as with a BH accretor. This estimate 
accounts for the specific lifetime of each binary in a po¬ 
tential EULX phase. We assumed that the typical den¬ 
sity of Milky Way Equivalent galaxies (MWEG) in the 















































4 



Fig. 2.— Upper limits on the number of EULXs within the 
sphere of radius R (solid for BH systems, dashed for NS systems) 
obtained under the assumption that a compact object can accrete 
at an arbitrary high rate (our reference model, Sec. |2.1[ |. Note 
that had we imposed the classical Eddington limit in our reference 
model, the predicted number of extreme ULXs would be zero. The 
dot-dashed and dotted lines correspond to the limited accr etion 
scenarios, with Lx calculated from the following models of |U07| 
and|S15] respectively. 


local universe is pmweg = 0.01 Mpc“^. Currently, the 
observations place the on ly confirmed EULX (H LX-1) at 
the distance of 95M pc (jWiersema et al.||2010[ but see 

dn 


Lasota et al.| (20151). Our predicted upper limit on the 
number ot EULXs at this distance is 90 and 93 binaries 
with BH and NS accretors, respectively. Both estimates 
significantly exceed the observed number of sources. 

Had we imposed the classical Eddington limit on the 
mass accretion rate in our reference model, the number 
of EULXs would be zero for both types of accretors. 

Table H introduces typical companion stars of our EU¬ 
LXs. Majority of the NS EULXs (92%) are found in 
binaries with Helium star donors: either a ^1.2 Mq 
H elium Hertzsprung gap star (HeHG) or ~1.8M0 He¬ 
lium Giant Branch (HeGB) star. The BH EULXs pre¬ 
dominantly contain a hydrogen-rich stars that have just 
evolved beyond main sequence. Typically, these are 
Hertzsprung gap (HG) stars (92.4%). The mass 
distributions of these most common EULX companions 
are presented in Fig. Table contains typical evolu¬ 
tionary routes that leM to the formation of EULXs^ _ 

As a result of implementing the slim disk model of|O07 


the X-ray luminosities of our systems drop by a factor o: 
at least 2.8. Nevertheless, in this framework it is still 
possible to significantly exceed the classical Eddington 
limit. The estimated number of the BH EULX systems 
decreases to 1 in 105 billion binaries, and the number of 
BH EULXs within 100 Mpc is 3 7 (Fig . [^. 

After employing the model of S15| the upper limit on 
a number of expected EULXs drops to 1 per 2 x 10^^ 
binaries. The expected upper limit on the number of 
EULX systems in the distance range of HLX-1 is 0.004. 
This is several orders of magnitude lower than in the two 
other models we tested. 

In all cases, the EULX luminosities are achieved during 




Fig. 3. — Mass distribution of the most common companions 
in NS (top) and BH fbottor n) sy stems during the EULX phase 
in our reference model (Sec. |2.1[ ). In the BH systems nearly all 
companions are on the Hertzsprung Gap (HG) with masses in the 
4-10 Mq range, while the NS companions typically enter the EULX 
phase as Helium Hertzsprung Gap (HeHG) stars and evolve to 
Helium Giant Branch stars (HeGB). 


the thermal timescale RLOF. The thermal RLQF mass 
transf er rate is calculated following [Kalogera &: Webbink| 
( 19961) as 


-^don 1 onLcion . . _i 

MRLOF,th =-= q X 10 . (5) 

where donor mass, radius and luminosity are expressed 
in solar units. For example an EULX with a BH can 
start with a ~ 10 Mq Hertzsprung gap donor, that has 
^ 20 Rq radius and luminosity of ^ 2 x IO^Lq. The 
donors in NS systems are evolved (post core He-burning) 
low-mass Helium stars that on the onset of RLOF have 
mass ^ 2Mq, radius ^ 30 Rq and luminosity ~ 2 x 
10 "^Lq. These parameters allow for very high RLOF mass 
transfer rates (MRLOF.th ~ 10“^-10 “^Mq yr“^. If this 
mass transfer is efficiently converted to X-ray luminosity 
(77 = 770), it will correspond to Lx ^ lO'^^-lO^^ ergs“^) 
alas on a very short timescale (rth lO^-lO'^yr). 

Typical evolutionary routes that lead to the formation 
of a potentia l BH and NS EULX when 77 = 770 are detailed 
in § |3.1H3.2[ The corresponding time behaviour of the 
mass accretion rate and maximum X-ray luminosity are 
shown in Figs. |^and[^ 

3.1. Black hole EULX 
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TABLE 1 

Compact accretors and their typical companions 


Accretor type'^ Companion type*^ 


NS EULX 

RG 

EAGB 

HeHG 

HeGB 

CO WD 

ONe WD 


2.6 X 10-3/MWEG 

6.35% 

6.4 X 10“3% 

48% 

44% 

0.5% 

0.8% 


BH EULX 

MS 

HG 

RG 

CHeB 

EAGB 

HeHG 

HeGB 

2.5 X 10-3/MWEG 

3.4% 

92.4% 

1% 

0.3% 

0.5% 

0.5% 

0.3% 


^ Type of accretor and the number of EULXs in the Milky Way Equivalent galaxy. 

Percentage of EULX systems with the same type of accretor, MS - Main Sequence, HG - Hertzsprung Gap, RG - Red Giant, CHeB - 
Core Helium Burning, EAGB - Early Asymptotic Giant Branch, HeHG - Helium Hertzsprung gap, HeGB - Helium Giant Branch, CO 
WD - Carbon-Oxygen White Dwarf, ONe WD - Oxygen-Neon White Dwarf. 


TABLE 2 

Typical EULX evolutionary routes 


Accretor / route 

%‘‘ 

Number/MWEG*’ 

Evolutionary route^ 

NS/l 

50% 

1.3 X 10“3 

MTl(2/3-l) MTl(8/9-l) SNl CE2(13-3/4;13-7) MT2(13-8/9) 

NS/2 

40% 

1.0 X 10-3 

MTl(2/3-l) SNl CE2(13-3/4;13-7) MT2(13-8/9) 

NS/3 

10% 

2.6 X lO-'^ 

Other 

BH/1 

97% 

2.5 X 10-3 

CE1(4-1;7-1) SNl MT2(14-l/2/3/4/5) 

BH/2 

3% 

6.8 X 10-5 

Other 


Percentage of the systems with the same type of accretor. 

Number of systems expected to be observed per Milky-Way Equivalent galaxy. 

Symbolic designation of the evolutionary routes; MT1/MT2 - mass transfer from the primary/secondary, SNl - supernova explosion. 
CE1/CE2 - common envelope phase started by the primary/secondary. Numbers inside the parentheses specif y the evolutionary phases 
of the stars: 1 - MS; 2 - HG; 3 - RG; 4 - CHeB; 5 - EAGB; 7 - HeMS; 8 - HeHG; 9 - HeGB; 13 - NS; 14 - BH (see|Belczynski et al. (|2008a|l 
for details). 


A typical BH EULX system evolves along the evo¬ 
lutionary route BH/1 presented in Table Its evolu¬ 
tion begins with a 33 Mq primary and 111^ secondary 
on an orbit with separation ~ 5500 Rq and eccentric¬ 
ity e = 0.56. In 5.5 Myr the primary starts crossing 
the Hertzsprung gap (HG). At that point the orbit has 
expanded to 5900 R© due to wind mass loss from the 
primary (now 30 M©). After about 10,000yr and signifi¬ 
cant radial expansion (1300 Rq), the primary begins core 
helium burning (CHeB). As the primary approaches its 
Roche lobe tidal interactions circularise the orbit. After 
some additional expansion (radius 1700 Rq) and mass 
loss the primary (18 Mq) initiates a common envelope 
(CE) phase. Following the envelope ejection the orbit 
contracts to 40 Rq and the primary becomes a Wolf- 
Rayet (WR) star with the mass of 11 Mq. After 6.2 Myr 
of evolution the primary undergoes a core collapse and 
forms a 7.2 Mq BH. At this time the orbit is rather com¬ 
pact (a = 47 Rq) and almost circular (e = 0.04). 

After the next 13 Myr the secondary enters the HG 
with a mass of 11 Mq and radius 10 Rq, and expands 
filling its Roche lobe (Riobe = 19 Rq). The mass trans¬ 
fer begins at orbital separation of 46 Rq. The luminos¬ 
ity of the donor is 21,000Lq. Initially, for a short pe¬ 
riod of time (6, OOOyr) the mass transfer proceeds on 
the donor thermal timescale with the RLOF mass trans¬ 
fer rate of Mrlof = 1-2 x 10“^ Mq yr“^ (corresponding 
to Lx = 5.8 X lO'^^ erg s“^). We allow the entire trans¬ 
ferred material to accrete onto the BH. After mass ra¬ 
tio reversal, the orbit begins to expand in response to 
the mass transfer and the RLOF slows down. However, 
for the next 2,OOOyr the donor RLOF rate stays above 


2.2 X 10“'* Mq yr“* (Lx = 10*^ ergs“*) and the system is 
still classified as a potential EULX. The evolution of the 
mass transfer rate throughout the RLOF phase is shown 
in Fig.ffl 

The RLOF terminates when the secondary has trans¬ 
ferred most of its H-rich envelope to the BH (with the 
final mass 14 Mq). The reminder of the secondary enve¬ 
lope is lost in a stellar wind and the secondary becomes a 
naked helium star with the mass of 2.5 Mq. At this point 
the binary separation is 230 Rq. The low mass helium 
secondary ends its evolution in the Type Ib/c supernova 
explosion forming a low mass NS (~ 1.1 Mq). The natal 
kick (if significant) is very likely to disrupt the system. 

3.2. Neutron star EULX 

A typical system (NS/I channel in Tablebegins as a 
10 Mq primary and a 5.6 Mq secondary on ^00 Rq orbit 
with an eccentricity e = 0.73. In 24 Myr the primary en¬ 
ters the HG. After 50, OOOyr of expansion on the HG the 
primary overfills its Roche lobe at the radius of 85 Rq . 
The orbit circularises and becomes 190Rq. We assume a 
non-conservative mass transfer in such a case, and allow 
50% of transferred material to accumulate on the main 
sequence secondary. At the end of this RLOF episode the 
primary becomes a low-mass helium star (2.2 Mq) and 
the secondary becomes a rejuvenated main sequence star 
(9.5 Mq). The separation has expanded to 640Rq. The 
low mass helium star evolves and expands after its core 
He burning is completed. When the radius of the primary 
(200Rq) exceeds its Roche lobe, the second phase of a 
non-conservative mass transfer is initiated. At this point 
the separation is 720 Rq and the masses are 2.06 Mq and 
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Fig. 4. — Typical evolution through Roche lobe overflow for a 
potential BH EUL X binary in our reference model (open circles, 
r] = rjo, see § |2.1[ |. A 11 MQHertzsprung gap star transfers its 
H-rich envelope to a 7Mq black hole. The mass transfer during 
the EULX phase is driven on a thermal timescale of the dono r 


The thick/red solid line and the hatched area represent, respec¬ 
tively, the mass accretion rate and the range of viewing angle de¬ 
pendent luminosities {8 = 0° —30° ) obtained for our Mrlof time 
evolution with the model of |S15| (/acc = 0.01, f'^ = 8). Simi¬ 
larly, the thick/red dashed line and thin/black dashed line show 
the m a.ss accretion rate and luminosity derived with the model of 
|O07 (/acc — 0.2-0.3, /■*" ~ 66). The classical Eddington limit is 
marked with dot-dashed line. 


9.45 for the primary and secondary, respectively. After 
a short phase of a RLOF (4, OOOyr) and after 28.6 Myr 
since the zero age main sequence (ZAMS), the primary 
explodes in electron capture supernova and forms a low 
mass NS (1.26 Mq). We assume zero natal kick in this 
case and the system survives the explosion. 

After additional 20 Myr the 9.3 Mq secondary leaves 
the main sequence and evolves through the HG to be¬ 
come a red giant. Expansion on the red giant branch 
leads to a RLOF that due to a large mass ratio and a 
deep convective envelope of the secondary turns into a 
CE phase. At the onset of the CE, the separation is 
600 R© and the secondary r adius is 330 Rf?^ . We assume 


energy balance for the CE ( 

Webbink 

1984), and obtain 

the post-GE separation of 50 Rq. J 

.’he secondary has 


lost its entire H-rich envelope and becomes a low-mass 
helium star ( 2 M 0 ). After core helium burning, the low 
mass secondary begins to expand and finally overfills its 
Roche lobe at the radius of 27 R© (the corresponding bi¬ 
nary separation is 55 R©). At this point the evolved he¬ 
lium star crosses the Helium Hertzsprung Gap and drives 
a high rate mass transfer (Mrlof = 7.8 x 10^ Mq yr“^; 
Lx = 8 X lO^^ergs”^) onto its NS companion. Since in 
our first approximation we have assumed a fully efficient 
accretion with no limit, the NS mass quickly increases 
to 2M0. The RLOF phase is very short (100yr), and it 
drives an extremely high mass transfer rate allowing this 
system to become a potential EULX with a NS accretor 
(see Fig. 1^. Finally, the secondary is depleted of its He- 
rich envelope and forms a naked CO core that cools off 
to become a CO WD. After 51.6 Myr since the ZAMS we 


Fig. 5. — Typical evolution through Roche lobe overflow for a 
potential NS EUL X bi nary system in our reference model (open cir¬ 
cles, 7] = 7]Q, see § |2.2| l. A 1.9 M© evolved helium star transfers its 
He-rich envelope to a neutron star. The mass transfer is driven on a 
very shor t thermal timescale (~ 100 yr) at a rate ~ 10~^ Mq yr~^ 
(see § |3.2[ |. The thick/red solid line and the hatched area represent, 
respectively, the mass accretion rate and range of viewing angle de¬ 
pendent luminosities obtained for our Mrlof evolution with the 
model of S15 (/acc = 0.01, ~ 8). Note that the model of |S15| 

considers tne case of BH accretion, however, in the case of a 
accretor, we assumed the same prescription for Lx- 

note the formation of a wide NS-WD binary (separation 
of 60 Rq), with the gravitational merger time exceeding 


the Hubble time {t„ 


= 2.8 X 101“^yr). 


4. DISCUSSION 

By allowing that all matter transferred through the 
RLOF is accreted {M^cc= ^klof) and converted effi¬ 
ciently into X-rays (Lx = where r] = rjo as in 

the standard disk), we are able to form a large number of 
potential BH and NS EULX systems. The EULX phases 
that we obtain are short (~ 10, OOOyr and ~ 100yr for 
BH and NS systems, respectively). Nevertheless, we ob¬ 
serve them in 0.1% of all simulated binaries. Our pa¬ 
rameter space covers all progenitors of X-ray binaries, so 
every 1 in 1, 000 XRBs should become an EULX during 
its evolution. Luminosities that we find in our simu¬ 
lations exceed lO'^^ergs”^. Such high luminosities are 
found both for BH and NS accretors. Particularly, the 
presence of the potential EULXs with NS accretors in 
our results seems to agree wit h the recent discover y of 
the NS ULX system M82 X-2 (jBachetti et al.||2014|. 

Evolution of ULX s with NS accretors w as also the topic 
of a recent work by Fragos et al. (2015). They used the 
BSE code for evolution of binaries and tire MESA code to 
calculate precisely the mass transfer phases. They found 
that NS ULX systems should exist in 13% of M82-like 
galaxies. They also found donors to be H-rich stars with 
masses in the 3-8 Mq range and 1-3 day orbital periods. 
The orbits in our NS EULX systems are wider (^ 30 days 
periods) and the companion stars are lighter (1-2 Mq) 
and of a diff erent type (ev o lved h elium stars) than those 
reported in Fragos et al. (2015). They obtained mass 
transfer rates up to A/rlof ~ 10“^ Mq yr“^ (see their 


log ( Mrlof / [Mq yr i] ) 
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Fig. 4), so they have reached mass transfer rates approx¬ 
imately as high as in our study. However, it needs to be 
noted that they have studied a much broader population 
of ULXs, while we have focused only on the brightest 
ones. 

A problem of the maximum X-ray luminosity available 
from a binary syst em was considered also by [Podsiad^ 
lowski et al. (2003). However, their study was confined 
to only 11 evolutionary routes, only BH binary systems, 
and quite limited parameter space (systems composed 
initially of a primary with mass in the 25 - 45 Mfr^ range 


and a 2 - 17 Mq secondary mass ran ge). Rappaport et al. 


( 2005 ) extended the simulations of Podsiadlowski et al. 


(2003) and were able to obtained ULXs with X-ray lumi- 
nosities of 3 x 10^^ ergs“^ (see their Fig. 11) for a 5 Mq 
BH and OMq donor and the case B mass transfer (HG 
donor). This is consistent with our typical BH EULX 
system. Their secondaries also hll their Roche lobe due 
to expansion of the envelope, and the ULX phase lasts at 
most a few Myr (about 10,000yr as EULX). However, 
their grid of parameters is far more sparse than ours, 
as they simulated only 52 specific evolutionary cases in 
c omparison to our 10^ evolu tion ary routes. 


Podsiadlowski et al. (2003) and Rappaport et al. (20051 
used detailed evolutionary code to obtained their results, 
while we have used simplified evolutionary formulae to 
cover larger parameter space. Both approaches have 
their advantages. With our approach we could not only 
confirm, but also extend the previously published results, 
both in terms of the accretor type (NS accretors possible) 
and in the range of potentially available mass transfer 
rates onto compact accretors in close binary systems. 

Reference model scenario leads to a significant overesti¬ 
mate of the number of the potential EULX systems in the 
local Universe. We find close to ^ 100 systems with a NS, 
and a similar number of systems with a BH within the 
100 Mpc radius. Note that these numbers should be con¬ 
sidered as upper limits. However, to date observations 
have re vealed only one sys tem with Lx > 10^^ergs“^, 
HLX-1 dFarrell et al.||20l'l ). 

Observational data indicates that HLX-1 is a transient 
system with recurrent outbursts on timescales of about 
1 year (e.g., Godet et al. 2009). Thermal-viscous in¬ 
stability mechanism was proposed to drive outburst in 
XRBs. However, this mechanism is not operational for 
mass transfer rates exceeding {a few) x i0“®MQyr“^ 
because the disk becomes too hot and thus constantly 
ionised (Lasota et al. 2015). If the periodicity observed in 
HLX-1 IS connected to thermal-viscous instability, then 
the mass transfer rate needs to be lower than the above 
threshold, and high X-ray luminosity is achieved by ef¬ 
fective beaming of radiation. 

Since our reference model provides only a very crude 
estimate of X-ray luminosity, we proceeded with inves¬ 
tigating the state-of-the-art global accr etion models of 
the super-Eddington accretion regime. 007 performed 
2-D radiation hydrodynamic simulations of supercriti¬ 
cal a ccreti on disks around BHs with the aP viscosity, 
while |S15| performed simulations of the magnetised ac- 
cretion disks in general relativity using radiation MHD 
code KORAL. The former simulations were fed by inflow¬ 
ing stream of gas circularising near i?cir = IOORgj while 
the latter were initialised as equilibrium torii threaded 


by seed magnetic field. Both models agree that there is 
signihcant mass loss in the accretion flow, driven either 
by radiation pressure itself, or by ra diatio n pressure and 
the centrifugal force. The work by |O07| provides a de¬ 
pendence of the mass accretion rate and luminosity on 
the mass input rate (here Mrlof)- However, this result 
depends strongly on the assumed location of the circu¬ 
larization radius, i.e., it implicitly assumes that no gas 
is lost outside R « 100i?G- In real RLOF systems, the 
outer edge of the disk is expected to be located much 
farther (up to about 2/3 R oche lobe radius of the accre¬ 
tor). The other approach (S15) was not limited by the 
disk truncation inside the simulation box, but rather by 
the computational time, which allowed flow to reach the 
inflow/outflow equilibrium only to radius i?eq ~ 100i?G- 
Inside this region, the gas flows out as wind down to 
I?wind ~ 20i?G, and shows roughly constant mass loss 
rate dMwind/di?. The total amount of gas lost in the 
system will therefore depend on the location of the outer 
disk edge, or rather the outer edge of the wind emitting 
region, through Eq. Because of poor understanding of 
the dynamics in the outer region of accretion disks, we 
assumed 7?wind/ R-out — 0.01. 

Within all of the discussed here accretion models bi¬ 
nary systems are able to breach the 10"^^ergs~^ EULX 
limit. Although, in the case of |S15| the probability of 
forming a binary EULX is a few orders of magnitude 
smaller than for the other cases. Limitation of the mass 
accr etion rate onto the BH in models developed in |O07| 
and |S15| results not only in lower X-ray luminosity, but 
also in different orbit evolution as compared to our typ¬ 
ical EULX case presented in Fig. Both effects (lower 
X-ray luminosity and different orbrt evolution) decrease 
the p redi cted number of EULXs in the local Universe for 
10071 and IS15I models. 

If disk winds are effective, matter which cannot be ac¬ 
creted is ejected from the system and takes away angular 
momentum. As a result the binary separation decreases 
and this prevents longer phases of high mass transfer. 
For the case of the |O07| outflow model, the total time 
spent in the EULX regime is about 3,000yr, as opposed 
to 10,000 yr in our reference scenario. As a result, we es¬ 
timated upper limit on th e number of the potential BH 
EULXs drops in the |O07| model from ^ 100 to 37 within 
the 100 Mpc radius. In Fig. |^we show how the Mrlof 
of our typical BH EULX (the reference model ) tran slates 
into Lx and Mace derived from the model of JO07 

When we limit the Mace according to Eq. |^S15), we 
find considerably fewer EULXs within the 100 Mpc ra¬ 
dius than in the case of our reference model (factor of 
^ 5 X 10^ fe wer). For accretion limited case employing 
the model of 007 we obtain only a few times fewer num¬ 
ber of EULXs (factor of ~ 3 fewer). Particularly, in the 
case of our typical BH EULX system (Fig. we find 2 
orders of magnitude smaller luminosity in the beam than 
in the reference model. In Fig. we show the il^cc and 
a range of L y computed for 9 = 0 ° - 30° (Eqs. [^andffl 
respectively. S§,dowski et al.|[^15[ ) corresponding to the 
mass transfer history Mrlof of our typical BH EULX 
(the ref erenc e m odel) . 

Both 10071 and IS15I considered accretion disks around 
BHs. The case of a NS accretion is far more complicated. 























































A strong magnetic fields associated with NSs may disrupt 
the disk far away from the accretor. Even if the magnetic 
field is weak, the emission will be weaker as the accretor 
mass is lower. On the other hand, the NS accretion effi¬ 
ciency ( 77 ) is higher than that of BHs as the matter and 
photons do not fall under the event horizon. To date, no 
comprehens ive simulations o f such a case have not been 
performed. Ohsuga (2007b) performed 2D simulations 
for a very limited range ot initial parameters. They ob¬ 
tained supercritical accretion rates for NS accretors and 
provided information that the mass accretion rate is 20- 
30% of that on to the BH for the same mass-input rate. 
Thus, we have assumed NS Eddington limited accretion 
case in the Oshuga et al. models, but we have allowed 
super-Eddington accretion for S^dowski et al. models. 

We note that ha d we impose the logarh hmic scaling of 
Lx with Mace (e.g., Poutanen et al. 2007), we would have 
not obtained any systems with luminosities in the EULX 
regime. There exist factors that can improve this situ- 


nebulae observed around some of ULXs (Pakull & Miri- 

oni|2002 

Russell et al.|2011 

Moon et ai.|2Ull|) may sug- 


emission. The relation between Mace and Lx has not 
been derived from first principles, and advanced numer- 
ical models do not seem to confirm this relation (e.g., 
Ohg]jggJ|200(^ S^owgH et al ||2015 )^ 


Gladstone et al. "(2013) and Heida et al. (2014) pre- 
sented programs to search for the companion stars of the 
ULX systems in the optical and infrared bands, respec¬ 
tively. The former group investigated close (~ 5Mpc) 
ULXs and discovered 13 ± 5 optical counterparts among 
33 ULXs. The masses of the companions have large 
uncertainties and generally only the upper limit is pro¬ 
vided in the mass range 5.7-16.1 Mq. In a few cases 
the lower limit is also present in the mass range 8.3- 
14.7 MfT). These constrai ns are in agreement with our re¬ 
sults. Heida et al. (2014) investigated 62 close ULXs and 
discovered 17 potential counterpart candidates. Based 
on the absolute magnitudes most of them (11) are Red 
Supergiants. According to our results the most common 
companions of the ULX systems able to reach the EULX 
regime are 6 Mq Hertzsprung gap stars with BH accre¬ 
tors and 1-2 Mq low mass Helium stars for NS accretors. 
However, it is possible that the Red Supergiants that we 
found to be in minority among the EULX companions, 
are more typical in the case of standard ULX systems. 


Detailed investigation of the properties of the entire ULX 
population will be presented in our forthcoming paper. 


5. CONCLUSIONS 

We conducted a proof-of-concept study to investigate if 
a binary system can form an ULX with the extreme mass 
transfer rate potentially able to lead to the X-ray lumi¬ 
nosities in the > 10"^^ ergs“^ range. This is at least hun¬ 
dred times more than what is exp ected for the Eddington 
limite d stellar-origin 100 Mq BH (Belczynski et al.|2010 
2014). Observations of HLX-1 in ESU 243-49 with X-ray 
luminosity of 1.1 x lO"*^ ergs“^ encouraged us to look into 
the problem. 

We find several evolutionary channels that lead to 
phases of a very high mass transfer rate in close RLOF 
binaries. These evolutionary phases can be extremely 
short, but it appears that many binaries experience such 
phases. The mass transfer rate may so high that (if Ed¬ 
dington limit, or any similar limit is breached) the X- 
ray luminosity may reach well above 10^^ergs“^. We 
have adopted two physical accretion disk models in X- 
ray binaries, and we have shown that for each model 
(that takes into account mass loss from accretion flow 
and photon trapping) EULX sources are formed within 
typical stellar populations. We note that increasing the 
geometrical beaming of radiation, many more binaries 
than considered in our study could possibly become ULX 
or EULX. 

It is found that about half of this potential extreme 
ULX systems host not BH, but NS accretors. This con¬ 
tradict the present consensus in the community, as the 
brightest ULX systems are commonly suspected to host 
IMBH primaries. Or at the minimum with stellar-origin 
BHs. 
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